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ABSTRACT 



1. INTRODUCTION 



The motions of the plasma and structures in and below 
the solar photosphere is not well understood. The results 
obtained using various methods cannot be in general con- 
sidered as consistent, especially in details. In this contri- 
bution we show a summary of the results obtained by the 
method we have developed recently. 

To study the photospheric dynamics we apply the local 
correlation tracking algorithm to the series of full-disc 
Dopplergrams obtained by Michelson Doppler Imager 
(MDI) aboard the SOHO satelite. The dominant structure 
recorded in Dopplergrams is the supergranulation. Under 
the assumtion that the supergranules are carried by the 
flow field of the larger scale, we study properties of this 
underlying velocity field. The methodology consists of 
an extensive data processing of primary data in order to 
suppress disturbing effects such as p-modes of solar oscil- 
lations or instrumental issues. Aditional coordinate trans- 
formations are also needed to make the data suitable for 
tracking. 

We perform comparative tests with synthetic data with 
known properties and with results of the time-distance 
helioseismology with a great success. Correlation coefi- 
cients of the comparison of mean components of the flow 
field are larger than 0.8, for the comparison of details in 
the vector velocity field the correlation coeficient is larger 
than 0.6. The results of the method applied to real data 
agree well with well-known features detected in the pho- 
tospheric velocity fields and reported by many other au- 
thors. A few case studies are shown to demonstrate the 
performance of the method. 

As a conclusion let's answer the question in the title. We 
believe that tracking of supergranules makes a perfect 
sense when studying the large-scale flows in the solar 
photosphere. The method we demonstrate is suitable to 
detect large-scale velocity field with effective resolution 
of 60" and random error of 15 m s"^ We believe that our 
method may provide a powerful tool for studies related 
to the dynamic behaviour of plasmas in the solar photo- 
sphere. 

Key words: Solar photosphere; velocity fields; local cor- 
relation tracking. 
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The solar photosphere is a very dynamic layer of the so- 
lar atmosphere. It is strongly influenced by the underly- 
ing convection zone. Despite years of intensive studies, 
the velocity fields in the solar photosphere remain not 
very well known. The evidence of the vigorous and some- 
times chaotic character of the motions of observed struc- 
tures in the photosphere (sunspots, granules, and other 
features) came already from the first systematic studies 
made in 19th century, of which let us at least mention 
the discovery of the solar differential rotation ( Camngton] 
|1859) . Motions in the photosphere are strongly coupled 
with magnetic fields. The large-scale velocity fields are 
very important for studies of the global solar dynamo. 

An attempt to describe the diff'erential rotation by a 
parabolic dependence did not make for clear results. Co- 
efficients of the parabola differed with traced objects and 
also changed with time, when tracing one type of object 
(reviewed e. g. by Schroterj 1985 1. The character of the 
differential rotation has never been in doubt. 

It follows from these arguments that a temporally variable 
streaming of the plasma exists on the surface of the Sun, 
which can be roughly described by the differential rota- 
tion. This streaming has a large-scale character, large- 
scale plasma motions were studied for example on the 



basis of tracking the magnetic structures (e. g. Ambroz 
[2001 a'bV The long-term Doppler measurements done by 



the MDI onboard SOHO make it possible to extend the 
studies of large-scale velocities in the solar photosphere. 
The knowledge of the behaviour of velocities in various 
periods of the solar activity cycle could contribute to un- 
derstanding the coupling between the velocity and mag- 
netic fields and of the solar dynamo function. 

There are at least three methods calculating photospheric 
velocities: 



1. Direct Doppler measurement - provides only one 
component (line-of-sight) of the velocity vector 
These velocities are generated by local photospheric 
structures, amplitudes of which are significantly 
greater than amplitudes of the large-scale velocities. 
The complex topology of such structures compli- 
cates an utilisation for our purpose. Analysing this 
component in different parts of the solar disc led to 
very important discoveries (e. g. supergranulation - 
lEtl[T954||Leighton et al.|[T9fel l. 



2. Tracer-type measurement - provides two compo- 
nents of the velocity vector. When tracing some pho- 
tospheric tracers, we can compute the local horizon- 
tal velocity vectors in the solar photosphere. Track- 
ing motions of sunspots across the solar disc led to 



the discovery of the differential rotation ( Carrington 
[18591 ). 

Local helioseismology - provides a full velocity vec- 
tor. Although local helioseismology (e. g. |Zhao| 
et al.' 2001) is a very promising method, it is still 



in progress and until now does not provide enough 
reliable results. 



Since the photosphere is a very thin layer (0.04 % of the 
solar radius), the large-scale photospheric velocity fields 
have to be almost horizontal. Then, the tracer-type mea- 
surement should be sufficient for mapping the behaviour 
of such velocities. In this field the local correlation track- 
ing (LCT) method is very useful. 

This method was originally designed for the removal 
of the seeing-induced distortions in image sequences 
( November 1986[ l and later used for mapping the motions 
of granules in the series of white-light images (November 
& Simon 1988 ). The method works on the principle of 



the best match of two frames that record the tracked struc- 
tures at two different instants. For each pixel in the first 
frame, a small correlation window is chosen and is com- 
pared with a somewhat displaced window of the same 
size in the second frame. A vector of displacement is then 
defined as a difference in the coordinates of the centres of 
both windows when the best match is found (cf. Fig. [T]). 
The velocity vector is calculated from this displacement 
and the time lag between two frames. 

The method needs a tracer - a significant structure 
recorded in both frames, the lifetime of which is much 
longer than the time lag between the correlated frames. 
We decided to use the supergranulation pattern in the full- 
disc Dopplergrams, acquired by the MDI onboard SOHO. 
We assume that supergranules are carried as objects by 
the large-scale velocity field. This velocity field is proba- 
bly located beneath the photosphere, so that the resulting 
velocities will describe the dynamics in both the photo- 
spheric and subphotospheric layer The existence of the 
supergranulation on almost the entire solar disc (in con- 
trast to magnetic structures) and its large temporal stabil- 
ity make the supergranulation an excellent tracer. 

The resulting velocities cover the whole solar disc. Hence 
the velocity field also describes a motion of the supergran- 
ulation in the areas occupied by active regions or by mag- 
netic field concentrations. This fact allows the results to 
be used to study the mutual motions of substructures like 
sunspots, magnetic field in active regions, background 
fields, and the quiet photosphere. 

The supergranules are structures with a strong convection 
coupling. The mean size of the supergranular cell is ap- 
prox. 30 Mm, with the size and shape that may dependent 



on the phase of the solar cycle. Supergranules are quite 
stable with a mean lifetime of approx. 20 hours. The 
internal velocity field in the supergranular cell is predom- 
inantly horizontal with the amplitude approx. 300 ms 
Due to the horizontality of the internal velocity field, the 
supergranules can be observed in Dopplergrams on the 
whole solar disc except for its centre. 

We do not propose that this method is capable of mea- 
suring velocities of order 1 ms"', but we do expect that 
the large-scale velocities will have magnitudes at least 
one order greater We also have to take the largest- 
scale velocities into account like the differential rotation 
or meridional circulation, which have velocities of at 
least 10 ms"'. If for example we take differential rota- 
tion described by the formula w[deg/day] - 13.064 - 
1.69 sin^ b — 2.35 sin"* b, we have to expect a velocity ap- 
prox. 190 ms"' mb - 60° in the Carrington coordinate 
system. The main goal of our study (and the proposed 
method of mapping the velocities should be proxy for this 
purpose) is to separate the superposed velocity field into 
the components and to investigate their physics. 



2. TUNNING THE METHOD 



A recent e xperience with applying this method to ob- 
served data Svanda et al. ](^005) has shown that for the 
proper setting of the parameters and for the tuning of 
the method, synthetic (model) data with known proper- 
ties are needed. The synthetic data for the analysis come 
from a simple simulation (SISOID code = 5/mulated 
Supergranulation as Observed In Dopplergrams) with the 
help of which we can reproduce the supergranulation pat- 
tern in full-disc Dopplergrams. 



The SISOID code is not based on physical principles tak- 
ing place in the origin and evolution of supergranulation, 
but instead on a reproduction of known parameters that 
describe the supergranulation. Individual synthetic super- 
granules are characterised as centrally symmetric features 
described by their position, lifetime (randomly selected 
according to the measured distribution function of the 
supergranular lifetime), maximal diameter (randomly ac- 
cording to its distribution function) and characteristic val- 
ues of their internal horizontal and vertical velocity com- 
ponents (randomly according to their distrib ution func- 



tion). For details about the SISOID code see ^Svanda et al. 



2006[ l. An example of the calculated synthetic full-disc 
Dopplergram can be seen in Fig. |2] 



3. METHOD OF DATA PROCESSING 



The MDI onboard SOHO acquired the full-disc Doppler- 
grams at a high cadence in certain periods of its opera- 
tion - one observation per minute. These campaigns were 
originally designed for studying the high-frequency oscil- 
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Figure 1. The sketch showing the performance of LCT method in detection of corresponding subframes in two frames 
separated by certain time-lag. 



lations. The primary data contain lots of disturbing ef- 
fects that have to be removed before ongoing processing: 
the rotation Hne-of-sight profile, p-modes of solar oscil- 
lations. We detected some instrumental effects connected 
to the data-tranfer errors. It is also known that the cali- 
bration of the MDI Dopplergrams is not optimal and has 
to be corrected to avoid systematic errors. While examin- 
ing long-term series of MDI Dopplergrams, we have met 
systematic errors connected to the retuning of the inter- 
ferometer. We should also take those geometrical effects 
into account (finite observing distance of the Sun, etc.) 
causing bias in velocity determination. In Fig. [3] the car- 
toon showing the individual steps in the data processing 
is provided. 

As input to the data processing we take a one-day observa- 
tion that contains 96 full-disc Dopplergrams in 15-minute 
sampling. Structures in these Dopplergrams are shifted 
with respect to each other by the rotation of the Sun and 
by the velocity field under study. 

First, the shift caused by the rotation has to be removed. 




For this reason, the whole data series (96 frames) is 
tracked with Carrington rotation rate, so that the helio- 
graphic longitude of the central meridian is equal in all 
frames and also equals the heliographic longitude of the 
central meridian of the central frame of the series. This 
data-processing step causes the central disc area ("blind 
spot" caused by prevailing horizontal velocity component 
in supergranules) in the derotated series to move with the 
Carrington rate. During the "derotation" the seasonal tilt 
of the rotation axis towards the observer (given by bo - 
heUographic latitude of the centre of the disc) is also re- 
moved, so that /7o = in all frames. 

Then the data series is transformed to the Sanson- 
Flamsteed coordinate system to remove the geometrical 
distortions caused by the projection of the sphere to the 
disc. Parallels in the Sanson-Flamsteed pseudocylindri- 
cal coordinate system are equispaced and projected at 
their true length, which makes it an equal area projection. 

The noise coming from the evolutionary changes in the 
shape of individual supergranules and the motion of the 
"blind spot" in the data series with the Carrington's rota- 
tional rate are suppressed by the k-co filter in the Fourier 
domain. The cut-off" velocity is set to 1 500 m s"' and has 
been chosen on the basis of empirical experience. 

The existence of the differential rotation complicates the 
tracking of the large-scale velocity field, because the am- 
plitudes and directions of velocities of the processed ve- 
locity field have a significant dispersion. We have found 
that, when the scatter of magnitudes is too large, veloc- 
ities of several hundred ms"' cannot be measured pre- 
cisely by the LCT algorithm where the displacement limit 
for correlation was set to detect velocities of several tens 
of ms Therefore the final velocities are computed in 
two steps. The first step provides a rough information 
about the average zonal flows using the diff'erential rota- 
tion curve 



oj - A + B sin^ b + C sin* b . 



(1) 



and calculating its coefficients. 



Figure 2. An example of the synthetic full-disc Doppler- 
gram produced by SISOID code. 



In the second step this average zonal flow is removed 
from the data series, so that during the "derotation" of 
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Figure 3. Cartoon showing the individual steps in the data processing. 



the whole series the differential rotation inferred in the 
first step and expressed by ([TJ is used instead of the Car- 
rington rotation. The scatter of the magnitudes of the mo- 
tions of supergranules in the data transformed this way is 
much smaller, and a more sensitive and precise tracking 
procedure can be used. 

The LCT method is used in both steps. In the first 
step, the probed range of velocity magnitudes is set to 
200 ms"', but the accuracy of the calculated velocities 
is roughly 40 m s ' . In the second step the range is only 
100 ms"' with much better accuracy. The lag between 
correlated frames equals in both cases 16 frame intervals 
(i. e. 4 hours in solar time), and the correlation window 
with FWHM 30 pixels equals 60" on the solar disc in the 
hnear scale. In one observational day, 80 pairs of velocity 
maps are calculated and averaged. 



4. METHOD CALIBRATION 



plotting the magnitudes of the model vectors versus the 
magnitudes of the calculated vectors (Fig. [4]i. The scat- 
ter plot contains more than 1 million points, and most of 
the points concentrate along a strong linear dependence, 
which is clearly visible. This dependence can be fitted 
by a straight line that can be used to derive the calibra- 
tion curve of the magnitude of calculated velocity vectors. 
The calibration curve is given by the formula 

Vcor = 1-13 Vcalc, (2) 

where Vcaic is the magnitude of velocities coming from the 
LCT, and Vcor the corrected magnitude. The directions of 
the vectors before and after the correction are the same. 
The uncertainty of the fit can be described by 1-cr-error 
15 ms"' for the velocity magnitudes under 100 m s"' and 
25 ms"' for velocity magnitudes greater than 100 ms"'. 
The uncertainties of approx. 15 ms ' have their main 
origin in the evolution of supergranules. 



In our tests we have used lots of variations of simple ax- 
isymmetric model flows (with a wide range of values of 
parameters describing the differential rotation and merid- 
ional circulation) with good success in reproducing the 
models. When comparing the resulting vectors of mo- 
tions with the model ones, we found a systematic offset 
in the zonal component equal to Voffset, zonal - -15 ms 
This constant offset appeared in all the tested model veloc- 
ity fields and comes from the numerical errors during the 
"derotation" of the whole time series. For the final testing, 
we used o ne of the velocity fie lds obtained in our previ- 

This field approximates 



ous work ( Svanda et al. 



2005 



the velocity distribution that we may expect to observe on 
the Sun. The model flows have structures with a typical 
size of 60", since they were obtained with the correlation 
window of this size. 

The calculated velocities (with Voffset, zonal - -15 ms~' 
corrected) were compared with the model velocities. Al- 
ready from the visual impression it becomes clear that 
most of vectors are reproduced very well in the direction, 
but the magnitudes of the vectors are not reproduced so 
well. Moreover, it seems that the magnitudes of vectors 
are underestimated. This observation is confirmed when 



^inodeled 



[m/s] 



Figure 4. Scatter plot for inference of the calibra- 
tion curve. Magnitudes of calculated velocities are 
slightly underestimated by LCT, but the linear behaviour 
is clearly visible. A line representing the 1:1 ratio is dis- 
played. The calibration affects only the magnitudes of the 
flows, while the directions do not need any correction. 
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5. COMPARISON TO TIME-DISTANCE HELIO- 
SEISMOLOGY 



effects of observations with different instruments or in dif- 
ferent conditions. 



Solar acoustic waves (p-modes) are excited in the up- 
per convection zone and travel between various points 
on the surface through the interior. The travel time of 
acoustic waves is affected by variations of the speed 
of sound along the propagation paths and also by mass 
flows. Time-distance helioseismology measurements iDu] 



vail et al. ( 1993 i and inversions provide (see references 



to Duvall's paper) a tool to study three-dimensional flow 
fields in the upper part of the solar convection zone with 
relatively high spatial resolution. 

The local correlation tracking method and the method of 
time-distance helioseismology provide surface or near- 
surface velocity vector fields. However, the results of 
these methods can be interpreted differently. While lo- 
cal helioseismology measures intrinsic plasma motions 
(through advection of acoustic waves), LCT measures ap- 
parent motions of some structures (granules or magnetic 
elements). It is known that some structures do not nec- 
essarily follow the flows of the plasma on the surface. 
For example, the supergranulation appears to rotate faster 

Some older stud- 



20001 



than the plasma ( Beck & Schou 
ies (see e. g. 'Rhodes et al. 1991 ' also suggest that the dif- 
ference in flow properties measured on the basis of struc- 
tures' motions and plasma motions is caused by a deeper 
anchor depth of these structures. An evolution of pattern 
may also play significant role (e. g. due to emergence of 
magnetic elements). 

Some attempts to compare the results of local helioseis- 
mology and the LCT method for large scales, with charac- 
teristic size 100 Mm and more, have been carried out by 
|Ambroz' (2005), but his results were inconclusive. The 
correlation coefficient describing the match of the ve- 
locity maps obtained by local helioseismology and the 
LCT method was close to zero. Nevertheless, there were 
compact and continuous regions of the characteristic size 
from 30 to 60 heliographic degrees with a good agree- 
ment between the two methods, so that one could not 
conclude that the results were completely different. In 
his study, many factors could be significant: the tech- 
niques were applied to different types of datasets (LCT 
was applied to low resolution magnetograms acquired at 
Wilcox Solar Observatory and the time-distance method 
used MDI Dopplergrams). Both techniques had very dif- 
ferent spatial resolution, and also the accuracy of the mea- 
surements was not well known. 

We decided to avoid these problems and analyze the 
same data set from the MDI instrument on SOHO. MDI 
provides approximately two months of continuous high- 
cadence (1 minute cadence) full-disc Dopplergrams each 
year. This Dynamics Program provides data suitable for 
helioseismic studies, and also for the local correlation 
tracking of supergranules. Thus, this is a perfect oppor- 
tunity to compare the performance and results of two dif- 
ferent techniques using the same set of data, and to avoid 



The selected dataset consists of 27 data-cubes from 
March 12th, 2001, 0:00 UT to April 6th, 2001, 0:31 UT, 
where each third day was used, and in these days three 
8.5-hour long data-cubes were processed (so that every 
third day in the described interval was fully covered by 
measurements). Each data-cube is composed of 512 
Dopplergrams with spatial resolution of 1.98" px^' in 
one-minute cadence. All the frames of each data-cube 
were tracked with rigid rate of 2.871 /irads"', remapped 
to the Postel's projection with resolution of 0. 12°px"' 
(1,500 kmpx"' at the center of the disc), and only a cen- 
tral meridian region was selected for the ongoing process- 
ing (with size of 256x924 px covering 30 heliographic de- 
grees in longitude and running from -54 ° to H-54 ° in lat- 
itude), so that effects of distortions due to the projection 
do not play a significant role. The time-distance inver- 
sion results were smoothed by a Gaussian with FWHM 
of 30 px to match the resolution to the LCT method, and 
only the horizontal components (v^, v,,) of the full veloc- 
ity vector were used. 

The results containing 27 horizontal flow fields were sta- 
tistically processed to obtain the cross-calibration curves 
for these methods. It is generally known that the LCT 
method slightly underestimates the velocities; thus, the 
results should be corrected by a certain factor. From the 
comparison of the jc-component of velocity we obtained 
parameters of a linear fit given by (numbers in parenthe- 
ses denote a cr-error of the regression coefficient) 



v.v, LCT = 0.895(0.008)v^, ,_d - 12.6(0.3) ms" 



(3) 



The correlation coefficient between v^^ lct and v^, t-d is 
p = 0.80. We assume that the time-distance measure- 
ments for V v_ t-d are correct and the magnitude of the LCT 
measurements, V j, lct, must be corrected according to the 
slope of Eq. ([3]). This correction factor has a value of 1.12, 
which is in a perfect agreement with the correction factor 
of 1 . 1 3 found in the tests of the same LCT code using syn- 
thetic Dopplergrams with the same resolution and similar 
LCT parameters. We assume that both velocity compo- 
nents obtained with the LCT method should be corrected 
by this factor The regression line of v, component is 



Vy, LCT = 0.56(0.0 l)v,, t-d + 0.4(0.2) m s" 



(4) 



After the slope correction using the Vx fits, the regression 
curve is slightly different: 



Vy, LCT = 0.63(0.01)V3, t-d + 0.4(0.2) ms" 



(5) 



with the correlation coefficient between v,., lct and v,,, t-d 
close to 0.47. The slope of the linear fit differs signifi- 
cantly from the expected value 1 .0. 

We have tested that this asymmetry is not related to the 
LCT technique. The tests did not show any preference in 
direction of flows measured by LCT or any dependence 
of the results on the size of the field of view (which does 



12* European Solar Physics Meeting, 8-12 Sep 2008, Freiburg, Germany 




50 100 150 200 250 50 100 150 200 250 

Figure 5. Left - Velocity field obtained by the LCT method. Right - Velocity field obtained by the time-distance technique. 
Both plots are centered at heliographic coordinates bo — 0.0°, /q = 214.3 °. Units on both axes are pixels in the data 
frame with resolution of 0.12 ° px^^ in the Postel's projection. 



not have a square shape in our case). In our synthetic data 
based test we did not encounter asymmetry of any kind. 

Instead of computation of the correlation coefficient of 
the arguments of both vector fields we decided to com- 
pute a magnitude-weighted cosine of A^. This quantity 
is given by 



Pw - 



2 Iv,-rfl 



|v,-rfl|vLC7-l ^ Z \yt-d\ cos Ay 

i:\y,-d\ SI ' ' 



(6) 



where \t-d is the time-distance vector field, v^cr is the 
LCT vector field and the summation is performed over 
all vectors in the field. The closer this quantity is to 1, the 
better is the agreement between both vector fields. Larger 
vectors are weighted more than smaller ones. We have 
found that in our case pw = 0.86. 

In addition to the detailed comparison of the vector fields, 
we compare the mean flows, the differential rotation and 
the meridional circulation. These flows can be quite sim- 
ply calculated from the results of both techniques. In both 
cases, they provide the mean zonal and mean meridional 
flows for the Carrington rotation No. 1974. The corre- 
lation coefficients are p = 0.98 for the zonal flow and 
p — 0.88 for the meridional flow. 



For a detailed comparison of the flow fields, we selected 
one data cube, representing 8.5-hour measurements cen- 
tered at 4:16 UT of Maixh 24th, 2001, /o=214.3 °. In this 
map, the correlation coefficient for the x-component of 
the velocity is p = 0.82, for the y-component p - 0.58, 
and for the vector magnitude: p = 0.73. The vector plots 
of the flow fields obtained by both techniques, shown in 
Fig. |5] in general seem to be quite similar to each other. 
However, many differences can be seen. The regions 
where the diff'erences are most significant correspond to 
relatively small (under 50 ms"') velocities. 



6. APPLICATION TO REAL DATA 



6.1. Long-term properties 

For the study of long-term evolution of surface flows, 
maps containing the mean zonal and meridional compo- 
nents were calculated. We can clearly detect that in the 
northern hemisphere the flow towards the northern pole 
dominates while in the southern hemisphere the flow to- 
wards southern pole prevails. The "zero line", the bound- 
ary between the flow polarity, is not located exactly on 
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the solar equator and seems to be shifted to the south in 
the period of increased solar activity (2001 and 2002). In 
agreement with Meunier et al. ( 1997| l and Cameron & 
[Hopkins ( 1998") we found the meridional flow stronger in 
the periods of increased solar activity by about 10 ms"' 
than in periods with lower magnetic activity. 

A similar map was made in the same way for the zonal 
component. The mean equatorial zonal velocity for all 
the data is 1900 m s"' . For all the processed data the de- 
pendence on latitude is close to a parabolic shape, the pa- 
rameters of which change slowly over time. The residua 
of the zonal velocity with respect to its parabolic fit given 
by 



(7) 



where b is the heliographic latitude and Vj, the mean zonal 
velocity in the given latitude, were calculated in order 
to see if we are able detect torsional oscillations in our 
measurements. The method reveals torsional oscillations 
as an excess of the mean zonal velocity with respect to 
the zonal velocity in the neighbourhood. The behaviour 
of torsional oscillations is in agreement with their usual 
description - the excess in magnitude is of the order of 
10 ms they start at the beginning of the solar cycle in 
high latitudes and propagate towards the equator with the 
progress of the 11 -year cycle. However, with our method 
the visibility of torsional oscillations decreases with in- 
creasing solar activity. In the periods of strong activity 
both belts are not so clearly visible since the large-scale 
velocity field and its parabolic fit are strongly influenced 
by the presence of magnetic regions. However, the tor- 
sional oscillation belts still remain visible when the mean 
zonal component is symmetrised with respect to the solar 
equator. We did not focus on the study of meridional flow 
or torsional oscillations depending on time and latitude, 
we just used them to check the ability and performance 
of our method. 

We investigated the coupling of equatorial zonal velocity 
(average equatorial solar rotation) and the solar activity 
in the near-equatorial area (belt of heliographic latitudes 
from -10 ° to -1-10 °). The average equatorial zonal veloc- 
ity incorporates the average supergranular network rota- 
tion and also the movement of degenerated supergranules 
influenced by a local magnetic field with respect to their 
non-magnetic vicinity. Indexes of the solar activity were 
extracted from the daily reports made by the Space Envi- 
ronment Center National Oceanic and Atmospheric Ad- 
ministration (SEC NOAA). Only the days when the mea- 
surements of horizontal flows are available were taken 
into account. As the index of the activity we have consid- 
ered the total area of sunspots in the near-equatorial belt 
and also their type. 

First, we computed the correlation coefficient p between 
the mean equatorial zonal velocity and the sunspot area in 
the near-equatorial belt and obtained a value of p = -0.17. 
We cannot conclude that there is a linear relation between 
these two indices. We find two different regimes which 
are divided by the velocity of approximately 1890 ms 



In one regime (77 % of the cases), the equatorial belt ro- 
tates about 60 ms"' faster (1910 + 9 ms"'; hereafter the 
"fast group") than the Carrington rotation, in the other one 
(23 %) the rotation rate is scattered around the Carrington 
rate (1860 ± 20 ms"'; hereafter the "scattered group"). 
The division in two groups using the speed criterion is ar- 
bitrary. If there are only two groups, they certainly over- 
lap and only a very detail study could resolve their mem- 
bers. The arguments for division in just two groups will 
follow. 

Detailed studies of the sunspot drawings obtained from 
the Patrol Service of the Ondfejov Observatory and the 
Mt. Wilson Observatory drawings archive revealed that 
in the "fast" group, the new or growing young active re- 
gions were present in the equatorial belt. On the con- 
trary, in the "scattered" group the decaying or recurrent 
active regions prevailed in the equatorial area. The de- 
celeration of the sunspot group with its evolution was no- 
ticed e.g. by ,Ruzdjak et al. (2004i. Moreover, our re- 
sults suggest that the new and rapidly growing sunspots 
in the studied sample (March to May 2001 and April to 
June 2002) move with the same velocity. This behaviour 
could be explained by emergence of the local magnetic 
field from a confined subphotospheric layer Accord- 
ing to the rough estimate (Corbard & Thompson 2002[ l 
the speed of 1910 + 9 ms corresponds to the layer of 
0.946 ± 0.008 /?0, where the angular velocity of rotation 
suddenly changes. During the evolution, the magnetic 
field is disrupted by the convective motions. 

The observed behaviour could be a manifestation of the 
disconnection of magnetic field lines from the base of 
the surface shear during the evolution of the growing 
sunspot group. This behaviour was studied theoretically 
by ISchtissler & Rempel| ( |2005p . They suggested a dy- 
namic disconnection of bipolar sunspot groups from their 
magnetic roots deep in the convection zone by upflow mo- 
tions within three days after the emergence of the new 
sunspot group. The motion of sunspots changes during 
those three days from "active" to "passive". The ac- 
tive mode is displayed by motions that are faster than 
those of a non-magnetic origin. The passive mode means 
mostly the deceleration of sunspot motion and influence 
of the sunspot motion only by the shallow surface plasma 
dynamics. The theory of the disconnection of sunspot 
groups from their magnetic root supports the division of 
the data set in two groups. 

We have also focused on how the presence of the mag- 
netic active areas will influence the average flow field. 
Since we found that the direct correlation is weak due 
to the existence of two different regimes, we decided to 
study the temporal change of both quantities. Our aim is 
to study whether an emerging active region in the near- 
equatorial belt will influence the average equatorial ro- 
tation. We computed numerical derivatives of the total 
sunspot area in the near-equatorial belt and of the aver- 
age zonal equatorial flow. We found that the correlation 
coefficient between both data series is p = 0.36 and is 
higher for the "fast group" (p = 0.41) than for the "scat- 
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Figure 6. Derivatives of the mean zonal velocity (solid curve) and the sunspot area in the near-equatorial region (dashed 
curve) in 2002. Both quantities correlate with each other well. 



tered group" (p = 0.24). The correlation is higher in pe- 
riods of increased magnetic activity in the equatorial belt. 
For example, for data in year 2001 the correlation coeffi- 
cient is P2001 = 0.58 and for year 2002 P2002 = 0.52 (see 
Fig. |6]l. In both cases, the correlation is higher for the 
"fast regime" (p ~ 0.7) than for the second group. 

Despite the apparent disagreement, our seemingly con- 
flicting results can be valid and are in agreement with 
results published earUer. As described by e. g. Hath- 



away & Wilson (1990 1 and explained in the model of 
Brun (2004|l, the solar rotation in the lower latitudes is 
slower in the presence of a magnetic field. In most cases, 
"spotty" equatorial belts seem to rotate slower than aver- 
age for the whole data series. However, it is clear that in 
most cases emerging active regions cause an increase in 
the rotation rate. This is in agreement with the generally 
accepted findin gs of Howard & Harvey ( |1970| l and Golub 



& Vaiana ( 1978f The relationship, obtained using a lin- 



ear fit on our data set, can be described by the equation 



Av ~ 0.2 AA„ 



sunspots 



(8) 



where Av is the change of the equatorial rotation speed 
with respect to the Carrington rotation and AA, 



sunspots 



the relative change of sunspot area in the equatorial belt 
(in 10"^ of solar hemisphere). We estimate that strong 
local magnetic areas rotate a few tens of m s"' faster than 
the non-magnetic surroundings. 



on October 7, 2004 at 16:30 UT at multiple wavelengths 
from ground and space instruments. 

At the point where the filament eruption begins (I - 56°, 
b = -26° in Carrington coordinates), we detected a steep- 
ening of the gradient in the differential rotation curve. 
During the eruption, the gradient flattens out and a dip 
forms. Although the differential rotation curves describe 
the mean zonal velocities on the full-disc, the change of 
its gradient signify the change in the stretch influencing 
the magnetic field in the loop over the area under study. 
The parameters of the smooth fitted curve did not change 
too much, however, the local residual with respect to the 
smooth curve changed at the latitude where the filament 
eruption starts. 

Fig. |7] displays the horizontal flows over a wide field 
of view measured using the our supergranule-tracking 
method and averaged over 3 hours, before and after the 
filament eruption. Before the eruption we can clearly see 
the north-south stream parallel to and about 10° east of 
the filament. This stream disturbs differential rotation and 
brings plasma and magnetic structures to the south. Al- 
though differential rotation tends to spread the magnetic 
lines to the east, the observed north-south stream tends to 
shear the magnetic lines. After the eruption, only a north- 
ern segment of the filament is visible and the north-south 
stream has disappeared. 



6.2. Flows around and beneath the erupting filament 

During three consecutive days of the JOP 178 campaign, 
Oct 6, 7, and 8, 2004, we observed the evolution of a fil- 
ament that was close to the central meridian. We also ob- 
served the photospheric flows directly below the filament 
and in its immediate area. The filament extends from -5° 
to -30° in latitude. A filament eruption was observed 



We measured the shear as a difference between the mean 
zonal component v\ in the area. We measured the av- 
erage zonal flow over boxes 2.3° on a side located 2.9° 
north and south of the point where the filament eruption 
appeared to start. Six 2-hour averaged of the flow fields 
were used to determine this figure. One can see that the 
shear velocity is increasing before the eruption and de- 
creasing after the eruption. One hour before the eruption 
the shear reached the value of (120±15) ms"' over a dis- 
tance of 5.2° (62 000 km in the photosphere). After the fil- 
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Figure 7. The 3-hour average of the flow field in the close vicinity of the starting point (I = 56°, b — —26°) before (left) 
and after (right) filament eruption. The filament observed by ISOON on October 7 2004 at 15.20 UT is superimposed. 
The dashed boxes denote areas used for the zonal shear calculation. 



ament eruption we observe ordinary differential rotation 
restoration below 30° south. 



6.3. Vertical structure of horizontal flows 

The issue of the our tracking-based method is that we 
measure displacements of supergranular structures in the 
series of processed Dopplergrams and interpret them as 
the large-scale velocity field in the solar photosphere. 
However, we cannot establish the range of depths in the 
solar convection zone, where the large-scale flows are 
well represented by the surface measurements. With the 
use of the data provided by the time-distance helioseis- 
mology we can determine this unknown parameter. We 
can also verify the assumption of our analysis: whether 
supergranules are subjected to the transport by the veloc- 
ity field on the larger scale. 

The large-scale flows are a combination of many types 
of motions of various spatial scales (such as differential 
rotation, meridional circulation, possibly giant cell con- 
vection, and others), which unfortunately cannot be reli- 
ably separated in components. These components may 
vary independently with time and depth. This implies 
that, in principle, when dealing with the large-scale flows, 
we cannot obtain unambigous results that do not allow 
any alternative interpretation. Nonetheless, the realistic 
numerical simulation should reveal the combined large- 
scale plasma flow, properties of which is comparable with 
the measured ones. Such numerical simulation, which 
is not present at the time, shall allow to distinguish var- 



ious components of the detected large-scale flow and to 
study them separately. The results therefore provide en- 
couraging large-scale flow properties, which are to be re- 
produced by models. 

For the comparison between the data obtained by our 
method and the results of the time-distance helioseismol- 
ogy we used high-cadence full-resolution full-disc Dopp- 
lergrams recorded in March and April 2001. In this pe- 
riod, 46 flow maps averaged over 8 hours of the vicinity 
of active region NOAA 9361/9393/9433 were processed 
by time-distance heUoseismology using the ray approxi- 
mation (see Zhao & Kosovichev 2004). From these maps 
we used six examples in six different days (March 3rd, 
March 28th, 29th, 30th, 31st, and April 25th). In these 
days, MDI data series did not contain many gaps, so the 
data were also suitable for our tracking method, and the 
field of view was far from the solar limb, so we can ex- 
clude any possible effect of undersampling in Doppler- 
grams near the solar limb. 

We had to slightly modify our method to match chal- 
lenges and properties of the time-distance data. Postel 
projection was used instead of Sanson-Flamsteed ones, 
because helioseismic data are computed in this geometri- 
cal projection (it conserves the main circles, therefore it 
is suitable for measurement of p-modes' travel-times). 

All the datasets were aligned with the centre of the field of 
view on Carrington coordinates / - 148.5° and b - 19°. 
The field of view has the total size of 512x512 pixels with 
the resolution of 0.92" px '. The datacubes from the lo- 
cal helioseismology contain 15 irregularly spaced depths 
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from 0.77 Mm to more than 80 Mm. The resuhs of our 
method have effective resolution of 60", but the results of 
time-distance helioseismology have 8". We smoothed he- 
lioseismic flow maps to match the resolution of 60". This 
also means that the structure of the internal flow in super- 
granulation is filtered out and does not disturb the study 
of the large-scale flow field, which is the most important 
deviation from the study of [Zhao & Kosovichev ( 2003^ . 
The velocity field obtained by our tracking method are in 
principle twodimensional - horizontal. The z-component 
of the time-distance flows suffers from effects of cross- 
talks, therefore we use only horizontal (i. e. in x and y) 
components of the time-distance flows. For each depth in 
the helioseismic datacube we calculated its similarity to 
the flow map obtained by tracking of supergranular struc- 
tures. As the measure of similarity we used the magni- 
tude weighted cosine of the direction difference. 

Our analysis shows that the large-scale flow does not vary 
much with depth down to 10 Mm in the quiet Sun regions 
(see Fig.[8]l. This fact means that the large-scale horizon- 
tal velocity field measured in the surface measurement 
represents well the large-scale horizontal dynamics in the 
layer of 0-10 Mm depth. The obtained results can be in- 
terpreted as a fulfilment of the basic assumption of our 
tracking method that supergranular structures may be in- 
deed treated as the objects carried by the large-scale ve- 
locity field and that the supergranules move in the depths 
down to 5 Mm coherently in the horizontal direction with 
a slight lost of coherence between 5 and 10 Mm. 

There is an alternative interpretation of the obtained re- 
sults. Assuming that the time-distance helioseismology 
is less reliable deeper down in the convection zone and 
that the large-scale flows do not vary much in the entire 
upper convection zone, our results put a level, where the 
reliability of the time-distance measurements is reason- 
able. In this interpretation, the time-distance flow mea- 
surements are reasonable in upper 10 Mm in the quiet 
Sun regions and upper 5 Mm in regions occupied by the 
magnetic field. There are recent papers (e. g. |Jackiewicz 
et al.!'2008) using another time-distance implementation. 



which suggest the range of reasonability of time-distance 
results only upper 3 Mm. 

To confirm our first interpretation, the structure of the 
supergranules itself must be analysed. Therefore we 
modified our analysis and instead of smearing the time- 
distance flow field to remove the signal of supergranu- 
lation internal mass flow, we subtracted the large-scale 
velocity field, so that only the internal velocities in su- 
pergranules remained in the flow map. From the maps 
containing only the separated supergranulation signal we 
calculated the horizontal divergence as the representation 
of down/up-flows. It is more stable than the measured 
z-component of velocity, which suffers of plenty from is- 
sues ( |Zhaol [2006) 1 ■ 

In the middle panel in Fig. [8] we see that the divergence 
signal is very similar in the quiet Sun region in depths 
0.77-6.4 Mm. Down to ~7 Mm the correlation is positive, 
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Figure 8. Upper: The similarity of the large-scale flows 
at various depths in the solar convection zone with the 
surface measurements by LCT. Middle: Correlation of 
the divergence signal representing supergranulation be- 
tween the small-scale flow maps at various dephts and 
the surface layers. Bottom: Trends of the opacity and 
superadiabaticity according to reference model S in the 
same level in the solar convection zone. 



implying the coherency in the structure of supergranules. 
Deeper down the correlation turns negative, which might 
suggest the evidence of return flow in supergranules. 

In the lowest panel in Fig. |8] the trend of the supera- 
diabaticity within the sub-surface layers of convection 
zone, coming from reference solar model S ( Christensen- 
Dalsgaard et al. 1996 1 is plotted. The superadiabaticity 
A* is defined as 



A* = 



1 d In p d In p 
7 d In r d In r 



(9) 



where y, p, and p are the state parametres (adiabatic ex- 
ponent, pressure, and density) of the plasma at distance 
r from the center of the Sun. Should the A* be negative, 
the layers at r are convectively unstable (e. g. discussion 
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in Demarque et al. 1997| l. It is to be noticed that A* turns 
more negative at the depths of 10-12 Mm, where it is 
assumed that the supergranulation should be formed, and 
remain very negative up to the solar photosphere, where it 
turns positive again (the photosphere is convectively sta- 
ble). It is the same layer, where the large-scale flow is 
coherent. It may be interpreted that the supergranules 
that form at depth some 10 Mm below the surface are 
carried by the large-scale velocity field, which operates 
here. Therefore the large-scale flow remain nearly con- 
stant within upper 10 Mm of the convection zone and is 
detected using our tracking method in surface measure- 
ments. 

In the same panel the trend of opacity, comming also 
from model S, is overplotted. Following the old opac- 
ity concept of conditions for onset of the convection 
( Schwarzschild 1975 1, we may conclude that from this 
point of view, there are suitable physical conditions at the 
depths of 10 Mm to form the convection pattern of super- 
granulation. 



7. CONCLUSIONS 

1. Tracking of supergranules makes a perfect sense 
when studying the large-scale field in the solar pho- 
tosphere and near subphotospherical layer. 

2. Our method provides results with effective resolu- 
tion of 60". 

3. Standard LCT lag between correlated frames is 
4hours, but is changeable, when averaging over 24 
hours, the noise level is at 15m/s (random error), 
other spans untested. 

4. SOHO/MDI provides enough material for some 
other studies, HMI/SDO will do better. 

5. We demonstrate the power of this method on a few 
case studies. We intend to continue in this research 
and expect the new results describing the amazing 
dynamics in the upper solar convection zone. 
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